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■ Abstract 

We examine the possibility of Bose-Einstein condensation (BEC) in two- 
dimensional (2D) system of interacting particles in a trap. We use a self- 
consistent mean-field theory of Bose particles interacting by a contact inter- 
action in the Popov and WKB approximations. The equations show that the 
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low T r the Bose-Einstein condensed state is not a consistent solution of the 
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equations in the thermodynamic limit. This result agrees with a theorem 



recently discussed by the author that shows that a BEC state is impossible 
for an interacting gas in a 2D trap in the thermodynamic limit. 
PACS numbers: 03.75.Fi, 05.30. Jp, 05.70.Fh, 32.80.Pj, 67.40.Db 
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^ | Recent experiments on alkali atoms [0] confined in three-dimensional (3D) 

hi . 

magnetic traps and cooled by evaporation techniques have led to the observa- 
tion of Bose-Einstein condensation (BEC). In these experiments the number 
of particles cooled has ranged from several thousand to millions. Experiments 
to date have naturally been in 3D systems. If experimentally it becomes pos- 
sible to make one dimension of the trap so narrow that the harmonic states 
are very greatly separated, then the system could be a reasonable simulation 
of a two-dimensional system. Further, there is the possibility of an adsorbed 
gas, such as spin-polarized hydrogen on liquid helium, M forming a 2D sys- 
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tern. 

The ideal Bose gas in two dimensions trapped in a harmonic potential has 
a Bose-Einstein condensation (BEC).|3], Q However, the author recently 
demonstrated how theorems by Hohenberg and Chester could be used to 
show that the trapped system can have no BEC, in the thermodynamic 
limit, if there are interactions that prevent the density from diverging at 
any position. The ideal gas has a divergence at the origin at and below the 
transition temperature, which excludes it from application of the theorem, 
and a BEC does occur. 

Our intent here is to test the general theorem by an explicit mean-field 
computation of the gas properties when interactions are present. We consider 
particles interacting by a contact potential. In 3D such an interaction is 
a pseudo-potential for a hard-core interaction. This is not the case in 2D, 
however it provides a simple model that can be analyzed by mean-field theory. 

Recently Griffin |j| discussed the generalized finite-temperature Hartree- 
Fock-Bogoliubov (HFB) equations of motion for a Bose gas in a nonuniform 
potential and the various approximations for treating it. He points out that 
a particularly useful gapless approximation is that of Popov. Giorgini et al 
(GPS) 0j have solved the equations in this form for three dimensions (3D) 
by use of the WKB approximation. Hutchinson et al |7j have solved the 
equations directly. 

We study the 2D system in the thermodynamic limit where the conden- 
sate equation may be treated exactly by the Fermi-Thomas approach. As 
in the ideal gas case, in the equations without condensate, there is a criti- 
cal temperature below which no value of the chemical potential satisfies the 
condition on the number of particles. At that point one usually invokes the 
presence of a condensate to satisfy the particle-number condition. However, 
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we find that there is no solution to the self- consistent equations for conden- 
sate and excited-state particles below the critical temperature. Apparently 
the noncondensed system becomes unstable at the critical temperature and 
makes a transition to some other state, but the state is not the BEC state. 
Possibly there is a Kosterlitz-Thouless transition, but we have not yet 
checked that hypothesis. [|| 

II. THE THERMODYNAMIC LIMIT FOR HARMONIC TRAPS 

Consider a 2D system of N Bose particles in a spherically symmetric 
harmonic trap. The harmonic potential is given by 

U(r) = \U^Y (1) 

where r 2 = x 2 + y 2 and R is a range parameter that will be handy for taking 
the thermodynamic limit. The angular frequency is 

where m is the particle mass. 

If we wish to take the thermodynamic limit we must increase the "volume" 
while keeping the average density fixed. The average density is proportional 
to p — N/R 2 where R is the range parameter in Eq. (jl|) chosen at some 
convenient temperature to be a distance within which the majority of parti- 
cles resides. Increasing the volume then implies weakening the potential, by 
increasing R, while N increases. From Eq. (^j we see that this requires keep- 
ing Nuj 2 = const while iV —>■ oo , This limiting process has been considered 
previously. [0, f|, [| Define a characteristic temperature 



where ks is the Boltzmann constant. We see that fcgTo = ^/NKlu remains 
constant as the thermodynamic limit is taken. For the ideal gas in a harmonic 
trap there is a phase transition^, |J at T c = T /y / ^(2) where ((a) is the 
Riemann ( -function. 

III. HARTREE-FOCK-BOGOLIUBOV EQUATIONS 

Recently Griffin || has discussed the derivation of a self-consistent mean- 
field treatment of the inhomogeneous interacting Bose gas valid at finite 
temperatures. One writes the field operator for the bosons as 

^ = $ + ^ (4) 

where $ = is the condensate wave function and t/> describes fluctuations. 
What results in the Popov approximation is a generalized Gross-Pitevskii 
equation valid at T > that now depends not only on the local condensate 
density n (r) = ($*$) but also on the density nr(r) = (ffi'ip) of particles 
in excited states. When we consider a contact interaction of strength g, the 
equation for the condensate is|| 

[A-<7n (r)]$ = (5) 

where the operator A is 

h 2 

A = -T^V 2 + U(r) -n + 2gn(r), (6) 

fi is the chemical potential, and n = no + n-r is the total local density. 

There is also an equation for ip that depends on uq and Ut- A Bogoliubov 
transformation of this equation leads to a pair of differential equations, which 
Hutchinson et al |7[] have solved in 3D by introducing an eigenfunction basis. 
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Giorgini et al M have used the WKB approximation to simplify and solve 
the equations for the excitations in 3D. The result of the same procedure in 
2D is the following: The excitation spectrum is 



e(p,r) = J A 2 - (gn y 



(7) 



where A = 2- + U(r) — /x + 2gn{v). The density of the excited particles is 



given by 
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» = J dp { [w 2 (p, r) + v 2 ( Pl r)] /(p, r) + w 2 (p, r)} 



with 



and 



/(P,r) = 
w 2 (p,r) 



v 2 (p,r) 



eP* - 1 
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The total number of particles satisfies 



N = J dv [n Q (r) + n T (r)] 



(8) 

(9) 

(10) 
(11) 

(12) 



The semiclassical approach makes sense if ksT ^> tkv. 

IV. A TRANSITION, BUT TO WHAT STATE? 

In 2D Eq. (^j) can be integrated. When one changes variables from p to 
y = f3e the integral takes on a particularly simple form: 

1 1 



n T lr 



dy 



1 



y 
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A 2 J^r=^~' a [ey - 1 
— {-In [l - exp(-v / t 2 ^ 2 )] + t - V^^ 2 } 



\/y 2 + s 2 



(13) 



where t = {3 (U(r) — /i + 2gn(r)), s = f3gno(r), and A 2 = h 2 j2 r nmkBT . 
When there is no condensate the density is simply 

4 >} W = -Tj In [1 - exp(-t)] (14) 
A 

We have solved Eq. fll4]) self-consistently for fi by numerical means and find 
that, just as in the ideal gas case, there is a solution only for T greater that 
some critical value. Sample results are shown in Fig. 1. Below the critical 
temperature we expect a new phase to exist. To see if the transition is to 
the BEC state we must consider the full set of equations. 



Baym and Pethick[10] have shown that the Fermi-Thomas approximation 
in which the kinetic energy is neglected is valid when N is large. In our 2D 
case one can show that the kinetic energy diminishes relatively as 1/N. Thus 
in the thermodynamic limit the Fermi-Thomas approximation is exact and 
Eq. (||) leads to 

n = -[//- U - 2gn T ] (15) 

This equation along with Eqs. ( |T2"D and (|TJ|) must be solved self-consistently. 
However, we see immediately that this is impossible because from Eq. ([15|) , 
s = t, and the exponential in Eq. (|l"3|) vanishes giving nonsense for all con- 
ditions. What has happened is that the lower limit in the integral form 
of Eq. (|H|) has vanished indicating that the long-wavelength phonons have 
destabilized long-range order just as in the homogeneous case. The BEC 
equations are inconsistent in agreement with the theorem discussed in Ref. 4, 
and there is no BEC in 2D in the thermodynamic limit. Note however that 
the mean-field equations for high temperature predict that there is a phase 
transition at some critical temperature. If those equations have any validity 
in describing a real system they tell us that the normal state of the interacting 
gas becomes unstable at some temperature. However, what state becomes 
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stable is not apparent from the present discussion, possibly a Kosterlitz- 
Thouless transition occurs. 

Real experiments are not done in the thermodynamic limit but with a 
finite number of particles. There can be a pseudo-condensation, a macro- 
scopic number of particles in the lowest state below some temperature that 
would go to zero in the thermodynamic limit. [|], [T]J Experiments in which 
this "transition temperature" is tracked as a function of N might be possible. 
Further theoretical computations for finite N are in order. 
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Figure Caption 

Fig. 1. a = —fi/ksT versus temperature divided by T for a sample set 
of parameters. Solid line: interacting gas with parameters 7(= gp 2 /kBTo) = 
1 and 7"o(= kBT p 2 m/h 2 ) = 1. Dotted line: Non-interacting case (7 = 
0, r = 1). With interactions there is a solution for \i above a critical reduced 
temperature, r c , but below, the self-consistent equations have no solution 
indicating there is no long-range condensate order. 
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